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Abstract 

In the present work we introduce a generalized lattice-gas model to study the multi- 
layer adsorption of interacting polyatomics on heterogeneous surfaces. Using an ap- 
proximation in the spirit of the well-known Brunauer-Emmet-Teller (BET) model, 
a new theoretical isotherm is obtained in one- and two-dimensional lattices and 
compared with Monte Carlo simulation. In addition, we use the BET approach to 
analyze these isotherms and to estimate the monolayer volume. In all cases, we found 
that the use of the BET equation leads to an underestimate of the true monolayer 
capacity. However, significant compensation effects were observed for heterogeneous 
surfaces and attractive lateral interactions. 
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1 Introduction 



The theoretical description of adsorption is a long-standing complex problem 
in surface science that presently does not have a general solution [Tf2"f3"P] . 

In 1918 Langmuir derived the monolayer adsorption isotherm kinetically for 
gas monoatomic molecules adsorbed on the homogeneous surface of adsor- 
bents without attractions among the adsorbed molecules [5]. Later, some 
theories have been proposed to describe multilayer adsorption in equilibrium 
PrTf8|9lfT0|llfl2fl3fl4fl5fl6fl7j . Among them, the one of Brunauer-Emmett- 
Teller (BET) [6] and the one of Frenkel-Halsey-Hill [7|8|9] are the simplest 
which provide the basis to construct more elaborate approaches. Those more 
elaborate analytic approaches take into account lateral interactions between 
the admolecules, differences between the energy of the first and upper lay- 
ers, surface energetic heterogeneity and so forth. These leading models have 
played an important role in the characterization of solid surfaces by means of 
gas adsorption. 

A more accurate description of multilayer adsorption should account for the 
fact that, in practice, most adsorbates are polyatomic. Even the simplest non- 
spherical molecules such as and O2 may adsorb on more than one site 
depending on the surface structure |17lll8|ll9f20ll21f22f23| . This effect, called 
multisite-occupancy adsorption, introduces a new complexity to the adsorp- 
tion theory. 

The difficulty in the analysis of the multisite statistics is mainly associated 
with three factors which differentiate the /c-mers statistics from the usual 
single-particle statistics. Namely, (i) no statistical equivalence exists between 
particles and vacancies; (ii) the occupation of a given lattice site ensures that 
at least one of its nearest-neighbor sites is also occupied; and (in) an isolated 
vacancy cannot serve to determine whether that site can ever become occupied. 

However, several attempts were done in the past in order to solve the k- 
mers problem. Among them, Onsager [23], Zimm [25] and Isihara [26] made 
important contributions to the understanding of the statistics of rigid rods in 
dilute solution. The FH theory, due independently to Flory [27] and to Huggins 
[28], has overcome the restriction to dilute solution by means of a lattice 
calculation. The approach is a direct generalization of the theory of binary 
liquids in two dimensions or polymer molecules diluted in a monomeric solvent. 
It is worth mentioning that, in the framework of the lattice-gas approach, the 
adsorption of /c-mers on homogeneous surfaces is an isomorphous problem to 
the binary solutions of polymer-monomeric solvent. 

More recently, two new theories to describe adsorption with multisite occu- 
pancy have been introduced. In the first, Ramirez- Pastor et al. [2"9~|3UII31] pre- 
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sented a model to study the adsorption of linear adsorbates on homogeneous 
surfaces. The model is based on exact forms of the thermodynamic functions 
of linear adsorbates in one dimension and its generalization to higher dimen- 
sions. In the second, which is called the fractional statistical theory of the 
adsorption of polyatomics (FSTA), the configuration of the molecule in the 
adsorbed state is incorporated as a model parameter [321133] . 

The treatments in Refs. [2^125 ,26,27,28,29.30.31 32 .331 are limited in their ap- 
plication because they are valid only for monolayer adsorption. There are few 
studies that take into account multisite occupancy at the multilayer regime 
|34f35f36f37f38f39f40] . Aranovich and Donohue [34"f35] presented a multilayer 
adsorption isotherm that should be capable to include multisite occupancy. 
Later, a closed exact solution for the multilayer adsorption isotherm of dimers 
was reported in Ref. [3D]. It demonstrates that the determinations of sur- 
face areas and adsorption energies from polyatomic adsorbate adsorption (also 
called /c-mers) may be misestimated, if this polyatomic character is not prop- 
erly incorporated in the thermodynamic functions from which experiments 
are interpreted. This treatment has recently been extended [37] to include the 
two-dimensional nature of the substrate. 

There are another two important physical facts which have not been suf- 
ficiently studied: 1) the effect of the lateral interactions between the ad- 
molecules and 2) the effect of surface heterogeneity in presence of multisite 
and multilayer adsorption. In the first recent paper [2H] extends the 

BET equation to include nearest-neighbor lateral interactions between the 
molecules adsorbed in the first layer. Following the configuration-counting pro- 
cedure of the Bragg- Williams approach and the quasi-chemical approximation, 
Ref. [38] gives us a simple statistical mechanical framework for studying mul- 
tilayer adsorption of interacting polyatomics. In the second case, multilayer 
adsorption of polyatomics on heterogeneous surfaces is a topic of interest in 
surface science since it comprises most features present in experimental situa- 
tions. Surfaces generally present inhomogeneities due to irregular arrangement 
of surface and bulk atoms, the presence of various chemical species, etc. In 
this context, a recent job shows how the monolayer volume predicted by BET 
equation differs from its real value when considering both the adsorbate size 
and the surface topography [3"D] . 

On the other hand, combined effects coming from lateral interactions, multisite 
occupancy and surface heterogeneity have been analyzed in the interesting 
paper by Nikitas [10] • In Ref. [ID], the author concludes that: (i) one can obtain 
an underestimation of the true monolayer capacity of the order of 25% when 
the adsorbate occupies more than one lattice site, and this underestimation 
will become worse if the effect of the multisite occupancy is coupled with 
heterogeneity effects and (ii) the attractive interactions in the gas adsorption 
lead always to a weak overestimation of the monolayer volume. 
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The previously discussed issues can be outlined in Fig. 1, where dimensions 
of complexity are sketched. In each axis we see a particular characteristic of 
the adsorption process that can (or cannot) be present in a particular model. 
Usually, the more features a model has, the more is its generality; and more 
complex the process it describes. In the multilayer regime, we can see models 
that include adsorption with lateral interactions, but they do not allow to 
multisite occupancy, or incorporate surface energetic heterogeneity, but they 
overlook lateral interactions in the adsorbed layer, etc. With Fig. 1 we try to 
show that, at present, some models venture in some of these dimensions of 
complexity but, the vast majority, lack the complete set of complexities. It 
is also clear that the farther away we are from the origin the more general 
(and usually the more complex) the model is. Even more, the upper plane 
(including multisite occupancy), is a set of models developed in the last decade 
that represent a leap forward in the lattice-gas theory. 

In this work we will attempt to take one step further away from the origin, 
into a more general and encompassing scheme, presenting a model that con- 
siders the totality of these dimensions of complexity (multilayer adsorption, 
lateral interactions, multisite occupancy, energetic heterogeneity and surface 
topography). For this purpose, a theoretical formalism is presented based upon 
the analytical expression of the adsorption isotherm with lateral interaction 
weighted by the characteristic length of the surface heterogeneous. In addi- 
tion, Monte Carlo (MC) simulations are performed in order to test the validity 
of the theoretical model. The new theoretical scheme allows us (1) to obtain 
an accurate approximation for multilayer adsorption on two-dimensional sub- 
strates accounting multisite occupancy, lateral interactions, energetic hetero- 
geneity and surface topography, and (2) to provide a simple model from which 
experiments may be reinterpreted. 

The present work is organized as follows. In Section II, we present the theo- 
retical formalism. Section III is devoted to describing the simulation scheme. 
The analysis of the results and discussion are given in Section IV. Finally, the 
conclusions are drawn in Section V. 



2 Model and theory 



In this section, we will present the model and make a revision of some existing 
theories describing the adsorption process in the framework of the lattice-gas 
approximation. We will focus on those contributions that give contextualized 
understanding of the work that we will present. 
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2.1 Model 

The substrate is modeled by a regular lattice of M sites with periodic boundary 
conditions, where the adsorption energy of the first layer e{ depends on each 
site i of the surface and the adsorbate is represented by fc-mers (linear particles 
that have k identical units). A A;-mer adsorbed occupies k sites of the lattice 
and can arrange in many configurations. This property is called multisite- 
occupancy adsorption. 

On the other hand, for higher layers, the adsorption of a k-mer is exactly onto 
an already absorbed one, with an adsorption energy of ke. Thus, the mono- 
layer structure reproduces in the remaining layers. The mechanism used to de- 
scribe the adsorption in the multilayer regime mimics the phenomenon called 
pseudomorphism. This was observed, by using low-energy electron diffraction 
technique [4"Tll42f43| . in the case of adsorption of straight chain saturated hy- 
drocarbon molecules on metallic surfaces. However, recent synchrotron X-ray 
scattering measurements show a different and more complex growth process 
in these systems [H]. Nevertheless, for simplicity, here we study a variant of 
the model that maintains the same mechanism to describe the adsorption in 
the higher layers. Because we analyze isotherms in the low coverage regime, 
we expect that this choice will not affect our results. 

Finally, attractive and repulsive lateral interactions are considered in the first 
layer and horizontal interactions are ignored in higher layers. In this sense, 
although the contribution from the secondary adsorption can already be sig- 
nificant, the density of the molecules in the second and higher adlayers is 
expected to be much lower than that in the first adsorbed layer [BET equa- 
tion can be applied at coverage not greatly exceeding (statistically) monolayer 
coverage] . Therefore, it seems to be satisfactorily enough to take into account 
only the interactions between the primarily adsorbed molecules [^5] . 

Under these considerations we can write the Hamiltonian as: 

M 



where the first term of the right-hand side (RHS) represents the adsorption 
energy of the Aq fc-mers adsorbed in the first layer (adjacent to de adsorbent) 
and the second term is the energy of the (N — Aq) /c-mers adsorbed on top 
of the first layer (second layer, third layer, and so on). The third and fourth 
terms correspond to the lateral interaction energy, where w is the interaction 
energy between two nearest-neighbor (NN) units belonging to different fc-mers 
adsorbed in the first layer (we use w > for repulsive and w < for attractive 
interactions); <7j is the occupation variable which can take the values if the 
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site % is empty or 1 if the site % is occupied and represents pairs of NN 
sites. Since the summation in the third term overestimates the total energy by 
including N%(k — 1) bonds belonging to the N% adsorbed k-mers, the fourth 
term subtracts this exceeding energy 

Our study will be restricted to the class of lattice-gas models in which the 
substrate is a regular array of individual adsorption sites where molecules 
can be deposited (also in a discrete manner). In the following sections we will 
review some analytical deductions of those models. A more in depth treatment 
is carried out in the original publications. 



2.2 Multilayer adsorption of non-interacting poly atomics on homogeneous 



One way to begin our analysis is to review the behavior of a system of identical 
particles, which will be adsorbed in a regular lattice of M identical sites of 
adsorption [H7]. The supposition is that, like in the BET model, adsorption 
is done on the surface of the solid or onto an already adsorbed particle. In 
this case, in stead of adsorbing spherically symmetric ad-atoms, we will use 
entities (that we will call k-mers) that occupy k consecutive lattice sites. 

Therefore, we will only have two possible unit-adsorption processes: (1) a k- 
mer occupies k consecutive empty surface sites; and (2) a k-mer adsorbs on 
top of an already adsorbed k-mer. This adsorptive process will form columns 
of k-mers on the solid. We must notice that a unit- adsorption where a k-mer 
adsorbs on top of two already adsorbed k-mers is prohibited. 

In the case of the unit-desorption processes, we will only be able to desorb a 
k-mer that is on the top of one k-mer column or a k-mer that is adsorbed on 
the solid surface, but has no other k-mers adsorbed on top of it. 

Under these conditions, the grand partition function of the system is: 



where n max (= M/k) is the maximum number of columns that can be formed, 
f2fc(ra, M) it is the total number of distinguishable configurations of n columns 
in M sites and £ is the grand partition function of a unique column of fc-mers 
that has at least one fc-mer in the first layer. 

On the other hand, we have that the grand partition function for the monolayer 
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(Hi) is: 



£ n k (n,M)X^ 



n=0 



(3) 



being n, in this case, the number of /c-mers adsorbed on the surface of the 
solid (first layer), and Ai the fugacity of the monolayer. VLk(n,M) is still the 
number of possible configurations with n k-mers in M sites. This quantity 
must be equal in Eqs. ()2]) and (J3]). 

If now we compare Eqs. ([2]) and ([3]) we can observe that they have a similar 
form. This allows us to write: 




M.T V^Ai J M T 



(4) 



If, in addition, inspired by this similarity, we propose the following ansatz: 
Ai = t- (5) 



Then 



Al = f = Z 1— => — = — 7T-, 6 

1 - p p Po 1 + c/Ai 



where c = q\/q = exp[— (5k{e^ —e)\ is the ratio between the partition function 
of a particle in the first layer and a particle in any other layer \fi = (/c^T)" 1 , 
being ks the Boltzmann constant]. And, then, we can write the monolayer 
coverage as: 




ains\ 

^ / M.T 



M.T 



Now, the total coverage (9) can be written in terms of the coverage of the 
monolayer (9\): 

1 - p/po 
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The theoretical procedure in Eqs. (jSj)-© provides the isotherm in the multi- 
layer regime from the isotherm in the monolayer regime. In fact: 

(1) By using 0\ as a parameter (0 < 61 < 1), the relative pressure is obtained 
by using Eq. OH]). This calculation requires the knowledge of an analytical 
expression for the monolayer adsorption isotherm. 

(2) The values of 0\ and p/po are introduced in Eq. (j9]) and the total coverage 
is obtained. The items (1) and (2) are summarized in the following scheme: 

0i + Ai(0i) + Eq. m^P/Po 

^9i+ p/po + Eq. © -> 9. 



Following the previous scheme, it is possible to obtain the exact multilayer 
isotherm for /c-mers in ID homogeneous surfaces. In fact, in Ref. [37] the 
expression for the monolayer coverage is proved to be: 



p/Po 



0i 



1 



(fc-i) 

k 



k-l 



kc(l - #i 



(fc-i) 

k 



9l 



k-l ' 



(10) 



Equations fl9]) and (|T0|) represent the exact solution of the ID model and, as 
it is expected, retrieve the BET equation for the case k — 1. 

Also, the previous scheme can be used to obtain an accurate approximation 
for multilayer adsorption on 2D substrates accounting multisite occupancy. 
In this case, the semi-empirical monolayer adsorption isotherm [33~P6] can be 
used 



di 



(fc-i)fli 



2(fc-l) /3 
(k Ul 



zkc{l-0 1 ) k + 1 



^0, 



(fe-i)fli 



(11) 



where ( is the connectivity of the lattice and z represents the number of 
available configurations (per lattice site) for a linear k-mer at zero coverage. 
Thus, z = 1 for k = 1 and z = (/2 for k > 2. 

Note that, for ( = 2, Eq. dTTJ is identical to the Eq. (HDD . Therefore, Eqs. © 
and ( II ip represent the general solution of the problem of multilayer adsorption 
in homogeneous surfaces with multisite occupancy. 

In the following section our aim will be to generalize this expression to include 
the lateral interactions following the methodology used in Ref. [38]. Clearly, 
the complexity of the isotherm is greatly increased, but the expression is still 
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a manageable one and in the case of zero interactions the previous case is 
retrieved. 



2.3 Multilayer adsorption of interacting polyatomics on homogeneous sur- 
faces 



Due to the complexity introduced in the analytical expressions because of the 
lateral interactions is that approximations are used to deal with this feature 
of the adsorption process. The most commonly used approximations are the 
mean-field approximation (MFA) [4"5f4"T] and the quasi- chemical approximation 

(qca) mm- 

As shown in previous work|47] . the configuration-counting procedure I" 2 ! of the 
QCA allows us to obtain an approximation that is significantly better than 
the MFA for polyatomics. Based on this finding, we restrict the rest of our 
discussion to the estimates obtained under QCA. 

In order to apply the theoretical scheme described in previous section, we start 
with the monolayer adsorption isotherm of interacting fc-mers adsorbed on a 
lattice of connectivity £ obtained from the formalism of QCA |47j . 



Ai 



9 1 exp ((3wz/2) 

kv(C,k) (f) 2(fc - 1] 



fc-i 



2k 



a 



z/2 



^ — (k — 1)6*1 



fc-i 



£(l-0i)-a 



fcC/2 



-.ft i 



where 



z = [2(C-l) + (A;-2)(C-2)] 
z( 9 1 (1-9 1 ) 



a 




2 



(13) 
(14) 

(15) 



and 



A = 1 - exp(-/3w). 



(16) 



2 The important assumption in this method is that pairs of NN sites are treated 
as if they were independent of each other (this assumption is, of course, not true, 
because the pairs overlap 
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Replacing Eq. (fl2|) into Eq. fl6]), we obtain 



1 + 



fc-1 r 



Id 



fcC/2 



Si exp (Pwz/2)(1 - [k-(k- 1)9^ 



fc-i 



2k 



a 



z/2 



,(17) 



Eqs. (Q and (IT7|) represent the solution describing the multilayer adsorption 
of interacting /c-mers on homogeneous surfaces in the framework of the QCA. 
This method is presented in more depth in [55] . 

Now we can see that we have moved upward in the plane of complexity and 
further more included interactions. The next section revises the role of surface 
heterogeneity and lateral interactions. 



2-4 Multilayer adsorption of interacting poly atomics on heterogeneous sur- 
faces 



In the two previous sections, we obtained the multilayer isotherm from the 
monolayer isotherm. It is possible demonstrate that this formalism still holds 
for interacting fc-mers and a given surface heterogeneity. However, this strat- 
egy leads to a complex solution that is not useful for practical purposes. To 
build a simpler function (easier to analyze), we have chosen to approximate 
the multilayer heterogeneous isotherm by a weighted sum of multilayer homo- 
geneous isotherms. 

The heterogeneous surface is modeled by two kinds of adsorption sites in the 
first layer (bivariate surface): strong sites with adsorption energy e{ and weak 
sites with adsorption energy e{- As seen later, these sites can be spatially dis- 
tributed in different ways (different topographies). Then, the total adsorption 
energy for an isolated /c-mer on the first layer with k\ monomers located over 
strong sites and k 2 monomers located over weak sites is 

Ei = k x e{ + k 2 e{. (18) 



Under these considerations, and using the formalism of the integral equation 
of the adsorption isotherm [TTj , the mean coverage 9 can be written as 

9 = ^2f(E i )9 loc (E i ), (19) 

Ei 



where f{Ei) is the fraction of fc-uples of ki strong sites and k% weak sites 
(ki + k 2 = k) and 6\ QC (Ei) represents the local multilayer adsorption isotherm 
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corresponding to an adsorptive energy E{. This local isotherm can be well 
approximated by using the multilayer adsorption isotherm associated to an 
homogeneous surface characterized by an effective value of c given by 

d = exp [-p (Ei - ke)] . (20) 



This value of q can also be expressed as function of c\ and C2, the values of c for 
homogeneous surfaces whose adsorption energies are e{ and e^, respectively. 
Thus, if the i-th term in Eq. (fT9l) corresponds to a k-mer with k\ units located 
over strong sites and ki units located over weak sites, then 




As an example, let us consider the multilayer adsorption of non-interacting 
dimers on a bivariate linear lattice, where strong and weak sites are spatially 
distributed in alternating homotattic patches of size I (I = 1,2,3, ■••). In 
this case, Eq. (fT9~j) has three different terms, being each one of them a dimer 
isotherm! 3 1 with a particular value of c, 




1 - 



1 -p/Po 
1 + (4ci - l)p/po 

1 - p/Po 



1/2' 



1/2" 



+ 



4:y/C^~ 1)P/P0 

1 - p/Po 
+ (4c 2 -l)p/Po 



1/2 



(22) 



The first [third] term in the RHS of Eq. $22 
a strong [weak] patch, on a pair of sites 



represents the adsorption within 



i) c i. 



with ci [c 2 ]. The 

fraction of (e[, e{) [{e{-, pairs on the lattice is (/ — 1)/2/ [ISj. The remaining 



[(4,4)L 



term of Eq. ( |22ll corresponds to a dimer isotherm on a (e{, £2) [or (£2) £ 'i)] 
pair (c = y/c\C2), being l/l the fraction of this kind of pairs on the lattice. 
As it can be observed, Eq. (122jl depends on I and the dimer isotherm sees the 
topography. 

In general, the number of terms in Eq. f lT9|) increases as the adsorbate size k 
is increased and this equation leads to a complex solution. This scheme can 
be notoriously simplified following the results in Ref. [39]. In this paper, the 
authors showed that: 



3 In the particular case of k = 1 and k = 2, it is possible to solve Eqs. ([9]) and (|10p 
to obtain single expressions of the multilayer isotherms corresponding to monomers 
and dimers. 
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If k ^> I (with k > 1), the multilayer adsorption isotherm can be represented 
by a single homogeneous isotherm 

= e Xoc (v^) • (23) 

For a topography where k <C I, the isotherm is 

= \o loc ( Cl ) + h loc (c 2 ). (24) 

The details of the topography are relevant only when k ~ I. In this case, it 
is possible to consider a simpler expression of the multilayer isotherm given 
by 



l -^f) #ioc fa) + (~) ^loc (v^) + V-yf) Oioc (c 2 ) . (25) 



Eq. (125]) (i) captures the extreme behaviors Eqs. ([231 and (1241) . and (22) 
approximates very well the complete Eq. (JT9l) in ID and 2D. In the 2D 
case, the local isotherm in Eq. (125]) is obtained from Eqs. fl9]) and ( flTT) . with 
C = 3, 4 and 6 for honeycomb, square and triangular lattices, respectively. 

Finally, we propose to extend Eq. f l25|) as to describe multilayer adsorption 
of interacting polyatomic molecules on heterogeneous surfaces. For this pur- 
pose, we propose to obtain the local isotherms from Eqs. ([9]) and ( IT71) . The 
advantages of using this simple description as a tool for interpreting multi- 
layer adsorption data and characterization of the adsorption potential will be 
shown in Section 4 by analyzing simulation results. 



3 Monte Carlo simulation 



The adsorption process is simulated through a grand canonical ensemble Monte 
Carlo (GCEMC) method. 

For a given value of the temperature T and chemical potential /i, an initial 
configuration with N /c-mers adsorbed at random positions (on kN sites) 
is generated. Then, an adsorption-desorption process is started, where each 
elementary step is attempted with a probability given by the Metropolis [49] 
rule: 

W = mm{l,exp[-P(AH - fiAN)]}, (26) 



where AH and AiV represent the difference between the Hamiltonians and the 
variation in the number of particles, respectively, when the system changes 
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from an initial state to a final state. In the process there are four elemen- 
tary ways to perform a change of the system state, namely, adsorbing one 
molecule onto the surface, desorbing one molecule from the surface, adsorbing 
one molecule in the bulk liquid phase and desorbing one molecule from the 
bulk liquid phase. In all cases, AiV = ±1. 

The algorithm to carry out one MC step (MCS), is the following : 

1) Set the value of the chemical potential \i and the temperature T. 

2) Set an initial state by adsorbing N molecules in the system. Each k-mer 
can adsorb in two different ways: i) on a linear array of (fc) empty sites on 
the surface or it) exactly onto an already adsorbed k-mer. 

3) Introduce an array, denoted as A, storing the coordinates of n e entities, 
being n e , 

n e = number of available adsorbed k— mers for desorption (n^) 



where n a is the sum of two terms: i) the number of fc-uples of empty sites 



on the surface and ii) the number of columns of adsorbed k-mera 4 [. 
4) Choose randomly one of the n e entities, and generate a random number 



4.1) if the selected entity is a fc-uple of empty sites on the surface then adsorb 
a fc-mer if £ < W™£, being W%£* the transition probability of adsorbing 
one molecule onto the surface. 

4.2) if the selected entity is a A;-uple of empty sites on the top of a column of 
height i, then adsorb a new k-mer in the % + 1 layer if £ < W^ l s k , being 
W^ l s k the transition probability of adsorbing one molecule in the bulk 
liquid phase. 

4.3) if the selected entity is a A;-mer on the surface then desorb the k-mer 
if £ < Wfag , being the transition probability of desorbing one 
molecule from the surface. 

4.4) if the selected entity is a k-mer on the top of a column then desorb the 
k-mer if £ < W^l, being W\^ lk the transition probability of desorbing 
one molecule from the bulk liquid phase. 

5) If an adsorption (desorption) is accepted in 4), then, the array A is updated. 

6) Repeat from step 4) M times. 

In the present case, the equilibrium state could be well reproduced after dis- 
carding the first m ~ 10 6 MCS. Then, averages were taken over ml ~ 10 6 MCS 
successive configurations. The total coverage was obtained as simple averages, 



Note that the top of each column is an available /c-uple for the adsorption of one 
k-mer. 



+ number of available k— uples for adsorption (n a ) 



(27) 




£e [0,1] 
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where (N) is the mean number of adsorbed particles, and (...) means the time 
average over the MC simulation runs. 



4 Results and discussion 

In the present section, we will analyze the main characteristics of the multilayer 
adsorption isotherms given by Eq. (1251) [with local isotherms obtained from 
Eqs. (jUJ) and (TTTj) ]. in comparison with simulation results for a lattice-gas of 
interacting k-mers on heterogeneous one- dimensional and square lattices. 

Heterogeneity is introduced by considering bivariate surfaces, i.e., surfaces 
composed by two kinds of sites in the first layer, strong and weak sites, with 
adsorptive energies e{ and e{,, respectively. Recent developments in the theory 
of adsorption on heterogeneous surfaces, like the supersite approach [5U], and 
experimental advances in the tayloring of nanostructured adsorbates [51)152] . 
encourage this kind of study. A special class of bivariate surfaces, with a 
chessboard structure, has been observed recently to occur in a natural system 
[53] . although it was already intensively used in modeling adsorption and 
surface diffusion phenomena [18 .54 .5 5155] . 

Bivariate surfaces may also mimic, to a rough approximation, more general 
heterogeneous adsorbates. Just to give a few examples, we may mention the 
surfaces with energetic topography arising from a continuous distribution of 
adsorptive energy with spatial correlations, like those described by the dual 
site-bond model [57], or that arising from a solid where a small amount of 
randomly distributed impurity (strongly adsorptive) atoms are added [5Sj. In 
both cases the energetic topography could be roughly represented by a random 
spatial distribution of irregular patches (with a characteristic size) of weak and 
strong sites. 

In the particular case studied in this article, the surface is modeled in two 
different ways: (1) as a chain of alternating patches of size / (see Fig. 2a); 
and (2) as a collection of finite homotattic patches in a chessboard-like array, 
where each patch is assumed to be a domain of equal size, I x I sites (see Fig. 
2b). In this model, the energy correlation length is simply given by the patch 
size. 

The computational simulations have been developed for one-dimensional chains 
of 10 4 sites, and square L x L lattices with L = 144, and periodic boundary 
conditions. With this lattice size we verified that finite-size effects are negligi- 
ble. Note, however, that the linear dimension L has to be properly chosen in 
such a way that it is a multiple of I. 
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We start analyzing what happens when the topography is changed. Fig. 3 
shows the behavior of the multilayer adsorption isotherms for k = 2, f3w = — 1 
and different topographies in ID as indicated. The energy difference between 
different patches has been chosen to be high (c\ = 1000 and C2 = 1) in order 
to emphasize the effects of the surface heterogeneity. We have identified the 
different topographies as lc for patches of size Z, and bp for the case of a surface 
with two big patches (/ —> oo). Symbols represent simulation results and lines 
correspond to theoretical data [Eq. (125]) ]. It can be seen that all curves are 
contained between the two limit ones: the one corresponding to lc and the 
one corresponding to bp. 

For 1 = 1, the adsorption energy of a dimer in the first layer is e[ + e{ for all 
configuration. In this condition, the system corresponds to a ID lattice-gas of 
interacting dimers on a homogeneous surface and, consequently, Eq. ( 123]) is 
exact. 

In general, for I > 1, this equation is approximate. For I = 2, the analytic 
isotherms agree very well with the simulation data. However, for p/po ranging 
between and 0.15, some differences between theoretical and numerical data 
are observed. This happens because Eq. (125]) has been built assuming that 
the three different pairs of sites are filled simultaneously and independently. 
However, for c\ 3> C2, the real process occurs in 3 stages: (i) the pairs of 
sites (s(,e{) are covered; (ii) the pairs (e^e^) begin to be filled and (Hi) the 
multilayer is formed. Note that in the first stage all the pair of sites (e{, e£) and 
(^2 5 ^i ) are removed. For this regime, a better approximation can be obtained 
by a semisum of two isotherms with c\ y c%. 

When I = 3, the agreement between the analytic isotherms and the simulation 
data is very good. In this case, for c\ 3> c<i the first stage does not eliminate 
all the pairs of sites (e{, ejjj) an d ( £ 2j £ i)j because each dimer occupies only 
two sites in the strong patches. For this reason, the range of validity of Eq. 
(I25p is wider than in the previous case. Now, if / = 4 or / = 5, the behaviors 
are similar to those observed for I — 2 or I = 3, respectively. In general, for 
even /, the first stage eliminates almost completely the pairs of sites (ef,^) 
and (^2, e{), while this does not happen for odd I. Finally, when / — > oo, 
the fraction of pair (e{, e^) and (e^^i) S oes to zero and Eq. (125]) is exact. 
This limit corresponds to the called large patches topography (bp surface in 
our model), where the surface is assumed to be a collection of homogeneous 
patches, large enough to neglect border effects between neighbor patches with 
different adsorption energies. 

We now analyze the effect of the lateral interactions on the behavior of the 
system. For this purpose, Fig. 4 shows the adsorption isotherms for k = 2, 
c\ = 1000, C2 = 1 and two different values of the lateral interactions (3w = — 1 
(attractive case) and (3w = 1 (repulsive case). In addition, for each value of 
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(3w, the limit topographies (lc and bp) have been considered (as seen in Fig. 
3, all curves corresponding to all topographies are contained between them). 

For repulsive couplings, the interactions do not favor the adsorption on the 
first layer and the isotherms shift to higher values of pressure. On the other 
hand, attractive lateral interactions facilitate the formation of the monolayer. 
Consequently, the isotherms shift to lower values of p/p and their slope in- 
creases as the ratio \(3w\ increases. In both cases (f3w = — 1 and f3w = 1), the 
agreement between theoretical and simulation data is excellent. 

From the curves in Fig. 4 (and from data not shown here for the sake of clarity) 
it is observed that: there exists a wide range of /3w's (—2 < f3w < 2), where 
the theory provides an excellent fitting of the simulation data. In addition, 
most of the experiments in surface science are carried out in this range of 
interaction energy. Then, the present theory not only represents a qualitative 
advance in the description of the multilayer adsorption of interacting fc-mers 
on heterogeneous surfaces, but also gives a framework and compact equations 
to consistently interpret thermodynamic multilayer adsorption experiments 
of polyatomics species such as alkanes, alkenes, and other hydrocarbons on 
regular surfaces. 

The effect of energetic heterogeneity (ratio between c\ and c 2 ) is analyzed 
in Fig. 5, where the degree of heterogeneity is varied by changing the value 
of c 2 between 1 and 100 with c\ fixed [c\ = 1000). As in Fig. 4, k = 2 
and f3w = — 1. Lines represent theoretical results and symbols correspond to 
simulation data (c 2 = 1: circles; c 2 = 10: triangles; and c 2 = 100: squares). 
For each set of values of the parameters, the limit cases corresponding to lc 
(open symbols) and bp (solid symbols) topographies are studied. As it can 
be observed from the simple inspection of the figure, the effect of topography 
is important in a range of c 2 /ci between 10~ 3 and 10 -2 and is practically 
negligible for 10 _1 < c 2 /ci < 1. 

To complete the discussion started in Fig. 3, we now evaluate the effect of the 
k-mei size on the adsorption isotherms. This study is shown in Fig. 6, where 
the lc and bp multilayer adsorption isotherms are plotted for (3w = — 1, 
Ci = 1000, c 2 = 1 and two different values of k {k = 2 and k = 10). One 
important conclusion can be drawn from the figure. Namely, the effects of 
topography and energetic heterogeneity tends to disappear as the size k is 
increased. 

The study in Figs. 3-6 was repeated for surfaces in 2D. The behavior of the 
curves (not shown here for brevity) is very similar to the one observed in ID. 

Summarizing, we have shown that just by using an expression of three terms, 
Eq. (1231) . we can approach very well the multilayer isotherm in ID and 2D for 
the adsorption of interacting polyatomics on heterogeneous surfaces. In the 
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next, we will use this approximation and Monte Carlo simulations to study 
how lateral interactions, multisite occupancy and surface heterogeneity affect 
the determination of monolayer volume predicted by the BET equation. 

In a typical experiment of adsorption, the adsorbed volume of the gas, v, is 
measured at different pressures and at a given fixed temperature, the total 
coverage is 9 = v/v m . Analyzing an isotherm with the BET equation 



it is possible to estimate the monolayer volume if we rewrite the previous 
equation as: 

-P/Po- (30) 



v(l-p/p ) cv m cu r 

This equation is a linear function of p/po- If we denote with a and b, the 
y-intercept and the slope of this straight line, respectively, we obtain 

(31) 



111 a + b 
and 



- + 1. (32) 
a 



The asterisk has been added in order to indicate that the quantities given by 
Eqs. ( 131]) and (132]) correspond to the prediction of the BET theory. Then, by 
means of a plot (the so-called BET plot) of the experimental data of v ^p° po ^ 
vs p/po, we can obtain an estimate of the monolayer volume and the param- 
eter c. Nevertheless, in the experiments it is commonly found that there are 
deviations from linearity in the BET plot. 

Following the scheme described in previous paragraphs, we carry out numer- 
ical experiments to determine, in different adsorption situations, how much 
the value of the monolayer volume predicted by the BET equation differs 
from its real value, v m . With this purpose, analytic and simulation isotherms 
were analyzed as experimental data. In this way, we have determined how ad- 
sorbate size, surface heterogeneity and lateral interactions, affect the standard 
determination of the monolayer volume. 

In Fig. 7 we plot the calculated value of v^/v m as a function of k for fixed 
values of c\ and C2 (namely c\ = 1000 and — 1) and different values of 
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the attractive lateral interaction. The repulsive interactions, not shown in the 
figure, show a marked deviation from the predicted BET equation, therefore 
wiping out compensation effects due to the surface heterogeneity. 

We can see three sets of plots in Fig. 7. Each set (1) is the result of the 
effect of different lateral interactions strengths: fiw = 0, circles; /3w = —0.5, 
squares and fiw = —1.0, triangles; and (2) depicts the limiting cases for the 
surface topography: open and solid symbols correspond to lc and bp surfaces, 
respectively. All other topographies must lay in between those two plots. 

If we analyze the behavior of the open symbols, we are analyzing, essentially, 
the homogenous case [39] where an already known feature for the homoge- 
neous scenario can be distinguished: the compensation effect of k is lower as k 
increases, as [37] stated for the non-interacting case. On the other hand, for bp 
surfaces, the compensation effect increases with k as seen in [35]. In all cases, 
for stronger (attractive) interaction strength the plots move upwards showing 
greater compensation effects. 

As shown in previous work [31)13?] for non-interacting /c-mers on homogeneous 
surfaces, the monolayer volume from the BET model diminishes with increas- 
ing values of the k-mer size. The data presented in Fig. 7 demonstrate that 
attractive lateral interactions and surface heterogeneity play a key role in the 
compensation of k-mer size effects. This finding is very important because 
most of the experiments in surface science are carried out in these conditions. 

Finally, the study in Fig. 7 is repeated for 2D surfaces (see Fig. 8). In this 
case, both the compensation effect for attractive interactions and the underes- 
timation of the monolayer volume for repulsive interactions (not shown here) 
are more important than in the ID case. The explanation is simple: in ID 
systems, particles will interact only at their ends, regardless of k; on the other 
hand, in 2D systems, each particle interacts at their ends, but also interacts 
along its k monomers. This makes the interaction energy grow linearly with 
k (if all first neighbors are occupied), and indicates that lateral interactions 
will play a more important role in two-dimensional adsorption systems than 
in one-dimensional ones. 



5 Conclusions 

In this work, we have studied the multilayer adsorption of interacting poly- 
atomic molecules onto heterogeneous surface. The polyatomic character of the 
adsorbate was modeled by a lattice gas of fc-mers. With respect to lateral in- 
teractions, the ad-ad couplings in the monolayer were explicitly considered in 
the solutions. The range of validity of analytical isotherms was analyzed by 



18 



comparing theoretical and MC simulation results. 

We also analyzed the ID and 2D BET plots obtained using the theoretical 
and simulation isotherms. For non-interacting fc-mers, we found that the use 
of BET equation leads to an underestimate of the true monolayer volume: 
this volume diminishes as k is increased. The situation is different for the case 
of interacting molecules over heterogeneous surface. Thus, attractive lateral 
interactions favor the formation of the monolayer and, consequently, compen- 
sate the effect of the multisite occupancy. In this case, the monolayer volume 
predicted by BET equation agrees very well with the corresponding true value. 
In the case of repulsive couplings, the lateral interactions impede the forma- 
tion of the monolayer and the BET predictions are bad (even worse than 
those obtained in the non-interacting case). Both the compensation effect for 
attractive interactions and the underestimation of the monolayer volume for 
repulsive interactions are more important for 2D systems. 
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Figure Captions 



Fig. 1: Schematic representation of the evolution of complexity in the theoret- 
ical adsorption models. 

Fig. 2: Schematic representation of heterogeneous bivariate square surfaces 
with chessboard topography. The black (white) symbols correspond to strong 
(weak) adsorption sites, (a) One- dimensional lattice and (b) square lattice. 
The patch size in this figure is / = 2. 

Fig. 3: Adsorption isotherms for dimers on one-dimensional lattices with /3w = 

— 1 and different topographies ID as indicated. The energy difference between 
different patches has been chosen to be high (ci = 1000 and c 2 = 1) in order 
to emphasize the effects of the surface heterogeneity. Solid lines and symbols 
represent theoretical and simulation results, respectively. For each set of values 
of the parameters, the limit cases corresponding to lc (open symbols) and bp 
(solid symbols) topographies are shown. 

Fig. 4: Adsorption isotherms for dimers on one-dimensional lattices with C\ = 
1000, C2 = 1 and two different values of the lateral interactions fiw = — 1 
(attractive case) and j3w = 1 (repulsive case). Solid lines and symbols represent 
theoretical and simulation results, respectively. For each set of values of the 
parameters, the limit cases corresponding to l c (open symbols) and bp (solid 
symbols) topographies are shown. 

Fig. 5: Adsorption isotherms for dimers on one- dimensional lattices with j3w = 

— 1. The degree of the surface heterogeneity is varied by changing the value 
of c 2 between 1 and 100 with c\ fixed [c\ = 1000). Lines represent theoretical 
results and symbols correspond to simulation data (c 2 = 1: circles; c 2 = 10: 
triangles; and c 2 = 100: squares). For each set of values of the parameters, 
the limit cases corresponding to lc (open symbols) and bp (solid symbols) 
topographies are shown. 

Fig. 6: Adsorption isotherms for dimers on one- dimensional lattices with f3w = 

— 1, Ci = 1000, c 2 = 1 and two different values of k (k = 2 and k = 10). Solid 
lines and symbols represent theoretical and simulation results, respectively. 
For each set of values of the parameters, the limit cases corresponding to lc 
(open symbols) and bp (solid symbols) topographies are shown. 

Fig. 7: Results of the BET plots for the adsorption in ID heterogeneous sur- 
faces with ci = 1000, c 2 = 1. Dependence on k of the fraction v^ a /v m for three 
different values of f3w: (3w = 0, circles; (3w = —0.5, squares and (3w = —1.0, 
triangles. For each set of values of the parameters, the limit cases correspond- 
ing to lc (open symbols) and bp (solid symbols) topographies are shown. 
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Fig. 8: As Fig. 7 for 2D surfaces: (3w = 0, circles; (3w = —0.1, squares and 
(3w = —0.5, triangles. Symbols connected by (solid) dotted lines correspond 
to results obtained from (theoretical) MC simulation isotherms. 
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